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Abstract 

We study the effects of switching social contacts as a strategy to control epidemic 
outbreaks. Connections between susceptible and infective individuals can be broken 
by either individual, and then reconnected to a randomly chosen member of the 
population. It is assumed that the reconnecting individual has no previous infor- 
mation on the epidemiological condition of the new contact. We show that recon- 
nection can completely suppress the disease, both by continuous and discontinuous 
transitions between the endemic and the infection-free states. For diseases with an 
asymptomatic phase, we analyze the conditions for the suppression of the disease, 
and show that -even when these conditions are not met- the increase of the en- 
demic infection level is usually rather small. We conclude that, within some simple 
epidemiological models, contact switching is a quite robust and effective control 
strategy. This suggests that it may also be an efficient method in more complex 
situations. 
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1 Introduction 



Mathematical models have been applied to the study of infectious diseases 
since more than a century ago. The last four decades have witnessed a burst 
of interest in quantitativ ely understanding the tra nsmission dynamics of a 
large number of diseases (lAnderson and Mayl . Il99ll ) . One of the key aims of 
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epidemiological mathematical models, and certainly the most relevant in terms 
of policy making, is the assessment of the effectiveness of control strategies to 
curb disease spreading. For many infectious diseases, the most widespread pre- 
vention measure is mass vaccination. However, if for a given disease vaccines 
are not known, or vaccination is not effective, other control measures have to 
be adopted. These mostly consist in some form of isolation of the infective in- 
dividuals. Isolation from an individual's social environment, however, can have 



seriou s psychological effects such as depression and stress (lHawryluck et al. 



20041 ). An even more drastic measure is quarantine, where individuals assumed 
to have been in contact with infective people are temporarily isolated from 
the rest of society. Daniel Defoe gives a forceful description of this practice 
and its tragic effects during the great plague of London (1665-1666) in his 
Journal of the Year of the Plague. A different strategy, that has helped pre- 
venting the spread of some sexually transmitted diseases, is contact tracing, 
where all the partners of an infective individual are located to be informed 
about the possibility of infection and, eventually, given adequate treatment. 
It has been argued, however, that this practice violates the right to privacy. 
In the ca se of HIV, for instance, it has sometimes met with considerable op- 
position (jBayer and Toomeyl . Il992l ). It must also be mentioned that all these 
measures require widesp read governmental action and, fr equently, allocation 
of substantial resources (lArmbruster and Brandeaul . 120071 ) . 



Recently, it has been suggested that a different kind of strategy, implemented 
by the individuals themse lves, could give surprisingly positive results in con- 
trolling disease spreading (IGross et al.l . l2006l ; IZanette and Risau-Gusmanl . l2008l ) 
The basic idea is that contact between non-infective and infective acquain- 
tances must be systematically avoided. This can be implemented by either 
non-infective or infective individuals, or by both. The effects of this simple 
strategy have been explored by means of network epidemiological models. In 
these models, each member of a population - hereafter, an agent- occupies a 
node of a network ([Keeling and Earned . 120051 ). Neighboring agents, which are 
connected by network links, can interact, and their epidemiological states (sus- 
ceptible, infective, recovered, etc.) thus change. For instance, the interaction 
of a susceptible agent and an infective agent may lead, with a certain prob- 
ability, to the infection of the former. The set of neighbors of a given agent 
represents the individual's acquaintances. 



In the framework of network models, the mechanism of avoiding contact with 
infective acquaintances amounts to breaking the links between susceptible 
and infective agents. If the social connectivity is nevertheless to be preserved, 
those broken links must be replaced by new connections. In the S IS (suscepti- 
ble — > infective — > susceptible) epidemiological model analyzed by lGross et al. 



(120061 ). it is susceptible agents who break contacts with their infective neigh- 
bors, what necessarily leads to the isolation of infective agents, and new con- 
nections are established only with other susceptible agents. Implicitly, this 
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assumes that the epidemiological state of all agents is known to any other 
agent, irrespectively of whether th ey are connected or not. More realistically, 



Zanette and Risau-Gusmanl (120081 ) consider links that are broken and recon- 



nected to randomly chosen members of the population. As described in the 
next section, we adopt here the same viewpoint. Furthermore, in order to 
avoid straightforward isolation of infective agents, in Section [3] we generalize 
the model in order to also allow infective individuals to switch connections 
from their susceptible neighbors to new contacts chosen at random. We show 
that, even under these less stringent assumptions on the mechanism of link 
switching, the outcome is still an effective reduction in the probability of infec- 
tion, which can eventually lead to the complete suppression of the disease. In 
Section HJ we study the case where the disease has an asymptomatic period, 
during which infective individuals ar e not perceived as such but appear as 
being susceptible, even to themselves (IFraser et all I2004J ) . In this situation, it 
may happen that an asymptomatic agent switches contact from an infective 
neighbor to a susceptible agent, thereby increasing the probability of transmis- 
sion. We find that, if the asymptomatic period is not very long, large enough 
switching rates can still suppress the infection. Results are summarized and 
discussed in the last section. 



2 SIS network model with switching links 



We consider a population formed by iV agents, situated at the nodes of a 
network. At a given time, each of them can be either susceptible (S) or infective 
(J). Initially, each agent has exactly k neighbors and, thus, the network has 
kN/2 links. The initial contacts are chosen at random, so that the resulting 
network is a random regular graph. The corresponding degree distribution, 
i.e. the distribution of the number of neighbors per agent, is thus a delta-like 
function centered at k neighbors. 

The infection is transmitted from infective agents to their susceptible neigh- 
bors, who become immediately infective (S — ► I), at rate A. Moreover, in- 
fective agents get cured at rate 7. Upon cure, they return to the susceptible 
state (/ — > S) and, subsequently, can be infected again. The diseases de- 
scribed by the SIS model are those that do not confer immunity against sub- 
sequent infections, and which have very short latent periods. Models of this 
class have proven particularly useful for the study of chlamydia and gonor 



Turner et all l2006h 



rhea (IHethcqte and Yorkd . Il984j ; iKretzschmar et all 119971 ; iGhani et al. 



1996 



Contact switching, in turn, is modeled as follows. All links joining susceptible 
agents (or, for short, susceptibles) with infective agents (or infectives) are 
broken at rate r, and the corresponding susceptible agent is then connected 
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to a randomly chosen agent in the population. Self-connections and multiple 
connections between pairs of agents are avoided. 



Taking into account the above dynamical rules, and using the pair approxima- 
tion discussed in the Appendix, it is possible to write down evolution equations 
for the variables n^(t) and rriAB(t), with {A, B} = {S, I}. They represent, re- 
spectively, the fraction of agents in state A and the fraction of links joining 
agents in states A and B. The equations read 



hj = -n/7 + Xm SI , 

fn SI = 2<ym n - <ym S i + XkKm SI —^- — - \m SI - rn s m SI , (1) 

ns 

ran = —27m// + Xm S i + 2XkKmsi — — . 

n s 

As the total number of agents and links is preserved, this system is completed 
by the relations mji + irtsi + mss = 1 and nj + ns = 1. The constant K = 
(k — l)/k derives from the pa ir approximation. The high-connectivity limit, 



K = 1 has been considered by 1 Gross et al.l (120061 ). 



In the remaining of this article we concentrate on the equilibrium properties of 
our system, obtained by equating the right-hand side of Eqs. (pQ) to zero. After 
some algebra, we find that the equilibrium fraction of infectives nj satisfies 
the fourth-degree polynomial equation 



m[-r(l - nif - (1 - nif - K 2 n] + 

+ K{1 -m){k\- (2 + X)m)} = 0. (2) 

The parameters f = r/7 and A = A/7 result from rescaling time in Eqs. flTJ), 
so that the time unit is the recovery time 7 -1 . The infection-free state, nj = 0, 
is a trivial solution of Eq. (j2J) for all values of the parameters. As expected, in 
fact, the infection cannot spontaneously appear in the absence of infectives. 
A second solution of the equation is always larger than unity. It therefore has 
no biological meaning, and will not be further discussed. As for the remaining 
solutions, it can be shown that for every value of the reconnection parameter 
f, there is a critical value of infectiveness A c at which a transition from the 
infection-free (A < A c ) to the endemic state (A > A c ), with persisting infection, 
occurs. 

To determine the nature of this transition we perform a stability analysis 
of the solutions of Eq. (j2J). When the reconnection parameter is below the 
threshold r c = [(2K — l)k + 1]/ (2k — 1), the situation is similar to the case 
without reconnection (f = 0). For small A the only stable non-negative solution 
is ni = 0. At A c , a transcritical bifurcation takes place. Namely, another real 
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solution becomes simultaneously positive and stable, whereas nj = becomes 
unstable. For larger values of the infectiveness, the equilibrium fraction of 
infectives grows continuously from zero. 



When, on the other hand, f > f c , the only real solution for A < A c is ni = 0. 
Now, the transition at the critical infectiveness A c is a tangent bifurcation, 
and two real and positive solutions appear simultaneously. For A > A c , the 
larger of these two solutions is stable and the other one is unstable, while 
the stability of nj = does not change. We have therefore a discontinuous 
transition from the infection-free state to the endemic state. Further increasing 
A, a transcritical bifurcation point is reached, where the unstable solution 
vanishes and exchanges stability with nj = 0. From then on, there is a single 
non-negative stable solution. Between the two bifurcations, the system has 
two stable equilibrium states, one endemic and the other infection-free. In 
this bistable region, the long-time asymptotic state is selected by the initial 
state of the population. An initial fraction of infectives larger than the unstable 
solution leads to the endemic state, whereas a smaller fraction leads to the 
extinction of the infection. The unstable solution thus plays the role of a 
threshold for the e pidemics to spread a nd, c onsequently, is sometimes called 
breakpoint density (I Anderson and Mayi . Il99ll ) . 



In Figure [T] analytical predictions of the equilibrium fraction of infectives, cal- 
culated from Eq. (T5]), are compared with numerical simulations of the epidemic 
process. As explained below, it makes sense to consider the solutions of Eq. 
(T5]) fixing the parameter K at both K = {k — l)/k (see the Appendix), and 
K — 1. For a reconnection parameter f = 0, the analytic al curve correspond- 
ing to K = (k — l)/k accurately fits the numerical results (jLevin and Durrettl . 
19961 ). When r ^ 0, on the other hand, for all but the smallest values of the 



reconnection parameter the curves that best fit the data are those with K = 1. 
A likely reason for this discrepancy lies in the fact that, depending on f, con- 
tact switching events generate networks with different degree distributions. It 
has previously been shown that, not unexpectedly, rewiring can lead to de- 



(Gross et al.. 


2006: 


Evans. 


2007) 



our model. It has also been pointed out that the pair approximation is not 
good for d escribing contact processes on networks with long-tailed degree dis- 
tributions (jPeyrard and Frand . l2005l ). The value K — (k — l)/k, in fact, results 



for the pair approximation applied to regular networks with k neighbors per 
node. A possible reason for the better fit with K = 1 could be that, with the 
progressive broadening of the degree distribution as rewiring proceeds, agents 
with many connections (k ^> 1) appear and become dominant in the epidemi- 
ological process. These agents, usually called super spreaders, are known to b e 



the main vector of propagation of many diseases (ILloyd- Smith et all 120051 ) . 
Taking K = 1 -which, as mentioned above, corresponds to the limit of k ^> 1 
for a regular network- is therefore a kind of effective approximation account- 
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Fig. 1. Equilibrium fraction of infective agents as a function of infectivity, for k = 3 
and four different values of the reconnection parameter f . Full and dashed lines are 
the analytical values of nj with K = (k — l)/k and K = 1, respectively. In the 
lower panels, the upper and lower branches of each curve respectively correspond 
to the stable and unstable equilibrium solutions. The lower branch represents the 
breakpoint density (see main text). Dots are results of single realizations in numer- 
ical simulations for a population of 500 agents (triangles) and 5000 agents (circles) , 
with 7 = 1. 

ing for the role of such agents. In the variants of the model considered below, 
again, better agreement between analytical and numerical results is found for 
K = 1. Consequently, all the analytical results presented in the following are 
obtained for this value of K. 

The possible asymptotic behaviors of our model can be summarized in a phase 
diagram on the parameter plane (A, f), as shown in Fig.[2]for three values of the 
average connectivity k, fixing K = 1. For each value of k, the plane is divided 
into three regions. In the leftmost region, corresponding to small infectivity, 
the fraction of infectives asymptotically vanishes and the population reaches 
the infection-free state. The central V-shaped region is the bistability zone, 
where the endemic state is stable but the infection-free state is still reached 
from initial conditions with sufficiently low levels of infectives. Note that, as 
discussed above, bistability does not occur for small reconnection parameters 
(f < f c ). Finally, for large infectivities, we have the region where the endemic 
state is the only stable equilibrium of the epidemiological process. To the left, 
this region is bounded by the transcritical bifurcation line, At = (f + l)/k. 
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Fig. 2. Phase diagram of the basic model, for three different values of k, and K = 1. 
For each value of k, the curve and the straight line represent, respectively, the 
tangent and the transcritical bifurcation. The V-shaped region between them is 
the bistability zone. The left-most region corresponds to the infection-free state, 
whereas the right-most region corresponds to the endemic state. 

It is interesting to notice that results similar to those derived from Eq. ([2]) can 
be obtained from the usual approach in mathematical epidemiology, where it 
is assumed that an epidemic outburst occurs as soon as an infective agent can 
generate more than one secondary infection (i.e., an infection event between 
two agents which involves a third, intermediary agent). This is quantified with 
the basic reproductive number Rq, which is the average number of secondary 
i nfections caused by a single infective agent in an otherwise healthy population 
(lAnderson and Mayl . 119911 ). The epidemic threshold is calculated by setting 
Rq = 1. In our case, a straightforward calculation of the basic reproductive 
number yields 



Rn 



kX 



r + A + 1 ' 



(3) 



from which the epidemic threshold turns out to be A^ = (f + l)/{k — 1). 
This is exactly the value at which the transcritical bifurcation occurs in our 
model when K = {k — l)/k. The slightly smaller result we obtain for K = 1, 
X T = (r + l)/k, is a consequence of the above discussed correction: it is well- 
known that wider degree distributions lead to smaller epidemic thresholds 
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( jAnderson and Mayl . Il99ll ; iPastor-Satorras and Vespignanil . l200ll ) . 



With the aim of making this basic model more realistic, two variants are 
considered in the following. In the first one, infective agents are allowed to 
reconnect their contacts with susceptibles. In the second, an asymptomatic 
stage in the disease is taken into account. 



3 Reconnection of infective agents 



Our basic model can be generalized to the more realistic situation where not 
only susceptibles are allowed to change their contacts, but also infectives can 
break their links with susceptible neighbours and redirect them to other mem- 
bers of the population. This possibility would stand for an altruistic attitude 
of infectives, in an attempt to inhibit the infection spread. To implement this 
variation, we introduce a new parameter p. When a link between a suscepti- 
ble and an infective agent is cut, with probability p it is the infective agent 
who keeps the connection, whose other end is redirected to a randomly chosen 
member of the population, avoiding self- and multiple connections. With the 
complementary probability, 1 — p, the connection is kept by the susceptible 
agent. Thus, this variation reduces to the basic model for p = 0. 

The evolution equations for the densities and the number of links -fixing, as 
advanced above, K = 1 in the pair approximation- are now 



hi = -tij7 + Ams/ , 
m 5 /=7(2m// - m SI ) + 

+Xkm S i-^^ — - Xm SI - rm SI [(l - p)n s + prij], (4) 
n s 

m II = -2 1 m II + Xm SI + 2Xkm SI ^ 

n s 

and the equation for the equilibrium fraction of infectives reads 



n/{-f[(l - n/) 3 +p(l - n/)(2n 7 - 1)] 

+ A(l-n/)(&-n/)-l} = 0. (5) 



As in the case of p — 0, this system undergoes a transcritical bifurcation 
where the infection-free state, with n 7 = 0, becomes unstable. The threshold 
infectivity at this bifurcation is = [f(l — p) + l]/k. The critical value of 
the reconnection parameter, above which a tangent bifurcation gives origin 
to two positive solutions, is r c = (k + 1)/(1 — p)(2k — 1). Note that, for f < 
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Fig. 3. Bifurcation diagrams for the fraction of infectives, as a function of the in- 
fectivity, for k = 10 and three different values of p, the probability that, for each 
susceptible-infective link, it is the infective who reconnects. From left to right the 
curves in each panel correspond to r = 10, 16 , 20 and 30. Full and dashed lines 
represent stable and unstable equilibrium solutions, respectively. 

f c , where the epidemic threshold is given by the transcritical bifurcation, \ T 
decreases with p. The reconnection of infectives, as expected, favors infection 
spreading. The disease is better controlled if only susceptibles are allowed to 
protect themselves against contagion. In the extreme case where only infectives 
reconnect, the infection threshold does not depend on f. 
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p 



Fig. 4. Bifurcation diagram for the fraction of infectives, as a function of p, for 
k = 10, A = 0.6 and several values of the reconnection parameter f. The upper 
and lower branches of each curve respectively correspond to the stable and unstable 
solutions. Dots are simulation results for a system of 5000 agents, with 7 = 1. 



For r > f c , the behavior of the epidemic threshold is more complicated. Figure 
[3] shows that the position of the tangent bifurcation is not monotonic with p. 
This leads to the unusual situation depicted in Fig. HI In certain regions of pa- 
rameter space, increasing p first makes the infection disappear, but a further 
increase of p leads to the reappearance of the endemic state through a second 
tangent bifurcation. For parameter values for which the infection is not sup- 
pressed, on the other hand, the fraction of infectives decreases monotonically. 



To check that these effects are not an artifact of the mathemati cal model, we 



have performed simulations of the corresponding contact process (jLevin and Durrett 



19961 ). whose results are represented by dots in Fig. HI We verify that simu- 
lations display the same qualitative behavior as the model, even though the 
analytical prediction is not as good as for p = (cf. Fig. [Q. This is proba- 
bly due to the fact that the degree distributions of susceptibles and infectives 
tend to be rather broad, as discussed above, and their respective parameters 
differ considerably. This difference grows with p, as infectives become, on the 
average, more connected than susceptibles. 

The phase diagram in the plane (A,f), shown in Fig. [5] for different values 
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Fig. 5. Phase diagram of the model with infectives who can reconnect, for different 
values of p. For each value of p, the curve and the straight line represent, respectively, 
the tangent and transcritical bifurcation. Compare with Fig. [2j 

of p, reveals the non-monotonic dependence on this parameter in the mutual 
crossings of the left boundary of the bistability regions. Note that, in any case, 
large enough reconnection rates always lead to the infection-free state. In fact, 
reconnection never increases the number of susceptible-infective links. In the 
next section, on the other hand, we consider a variation of the basic model 
where susceptible-infective links can be created by contact switching. 



4 Effect of asymptomatic agents 



Several common diseases possess an asymptomatic stage, intermediary be- 
tween the susceptible and the infective phases, during which individuals car- 
rying the infection do not exhibit symptoms, and are therefore seen as non- 
infective by themselves and by the rest of the population. In the asymptomatic 
stage, however, an infective person may already transmit the disease by con- 
tagion. It is well known that during the several plagues that ravaged Europe, 
infective people who did not yet have visible symptoms became one of the 
main vectors for the propagation of the disease, when they fled the cities. It 
is therefore not unreasonable to fear that, for a disease with an asymptomatic 
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phase, the strategy of contact switching might reinforce the spread of the 
disease, instead of curbing it. 

To address these issues, in this section we analyze a model with three types 
of agents, namely, S, A (asymptomatic), and /. It is basically an extension of 
the model studied in Section [21 because the interactions between susceptibles 
and infectives are the same. We assume that asymptomatic agents do not 
know that they are infective, and therefore act as susceptible individuals: they 
reconnect their links with infective neighbors at the same rate r as susceptibles. 
On the other hand, asymptomatic agents can infect their susceptible neighbors 
at the same rate A as infectives. Note that a main difference with the two 
models analyzed above is that, now, susceptible-asymptomatic links -which 
allow for the transmission of the infection- can be created by reconnection from 
asymptomatic-infective links, which join two agents carrying the infection. 

To characterize the transitions A —> I and / — > S we introduce a new param- 
eter a. The rate at which asymptomatic agents become (symptomatic) infec- 
tives is 7/ct, while the rate at which infectives become susceptible is 7/ (1 — a). 
With these definitions, the asymptomatic and infective phases last, on the av- 
erage, a/7 and (1 — a)/7 time units, respectively. The mean total infective 
period is, as before, 7 _1 . In the limit a — > 0, we recover the model of Section 
[2J whereas for a — > 1, the model is equivalent to the usual SIS model without 
reconnection. 

The evolution equations for this model are 



7 . 7 

rij = -- m + -n A , 

1 — a a 

7 

h A = n A + Xk(m SI + m SA )/2, 

a 

7 , n \ 1 7 ,fc rn SI + m SA 

m S i = - (2m 7 / - msi) + -m SA - X-m S i Xm S i 

1 — a a 2 ns 

-rm S i{n s + n A ), 

7 x , msi + m SA 

m ss = msi ~ Xmssk h rn s m S i, 

1 — a ns 

m n = -2 - 7 m u + -rriAi, (6) 
1 — a a 

7/ \ 7 . x k msi + m SA 

m A i = (2m AA - m AI ) - m AI + X-m S i h Xm S i 

a I — a 2 ns 

-rm AI (n s + n E ), 

7 ,,k m SI + m SA 

m AA = -2-m AA + X-m SA h Xm SA + rn A m AI . 

a 2 ns 

The system is completed by the equations n s + n A + n t = 1 and m S s + m SA + 
m S i + m AA + m AI + m n — 1. In spite of the rather imposing look of this math- 
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ematical problem, it is still possible to find a single (sixth-order) polynomial 
equation for the fraction of infectives nj at equilibrium. The equilibrium frac- 
tions of infective and asymptomatic agents are related by the simple equation 
(1 — a)riA = arii. Thus, at equilibrium, a exactly coincides with the fraction 
of agents carrying the infection which do not exhibit symptoms. This, in turn, 
makes it possible to find the equilibrium fraction of susceptibles, ns = 1— nj / a. 

As in the models studied in Sections [2] and [31 the infection-free state, nj = 0, 
is an equilibrium solution for any value of the parameters. The transcritical 
bifurcation, at which rij = changes its stability, occurs now at a critical 
infectivity At = [1 + (1 — a)r]/A;[l + (1 — a)ar]. As before, At grows monoton- 
ically with f. Here, however, this growth is not unbounded, and the critical 
infectivity tends to Aoo = (ka)' 1 for f — ► oo. In other words, for A > A^ 
no reconnection rate is large enough to suppress the disease. This can be un- 
derstood by inspecting the basic reproductive number corresponding to the 
present epidemiological process, 

iZo= - ^ + - kl . (7) 

A + a- 1 f + A+ (1 - a)- 1 



The two terms in the right-hand side of this equation arise from the number 
of secondary infections respectively generated during the asymptomatic and 
symptomatic phases. As expected, only the latter depends on the reconnection 
parameter. But, for A > A^, the first contribution is already larger than one. 
In other words, the infectiveness is so high that when symptoms appear - 
and reconnection is turned on- the infective agent has already caused enough 
secondary infections to trigger the spreading of the infection. 

A numerical analysis of the equilibrium solutions of the system in Eq. (J6]) 
reveals that, as before, a tangent bifurcation occurs at a critical infectivity 
smaller that Ay, for sufficiently large reconnection rates. Figure [6] shows the 
phase diagrams for three values of a. We see that the critical infectivity corre- 
sponding to the tangent bifurcation varies non-monotonically with f: it reaches 
a maximum A2 at an intermediate value of the reconnection rate, and then 
approaches a limit Ai as f — > 00, with Ai < A2 < At- 

Thus, the picture of the effects of reconnection is as follows. For large values of 
infectivity, A > Aoo, the endemic state is stable for all f. When the infectivity is 
in the interval (A2, At), increasing f takes the system from the endemic state to 
the bistable region, which here extends to arbitrarily large values of the recon- 
nection rate. Within this region, as discussed in Section [2J there is a threshold 
initial density of infectives necessary to sustain the epidemics. This breakpoint 
density grows with the reconnection rate and, for f — > 00, approaches a value 
smaller than the corresponding equilibrium fraction of infectives. Under these 
conditions, even an arbitrarily large reconnection rate is unable to suppress 
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Fig. 6. Phase diagram for a disease with an asymptomatic stage of duration a/7, and 
three values of a. For each of them, the left-most and right-most curves represent the 
tangent and transcritical bifurcations, respectively. Compare with Figs. [2] and [5j The 
arrows indicate the position of the vertical asymptotes of the curves for a = 0.2 
and a = 0.4: from right to left, they correspond to X\(a = 0.2), X±(a = 0.4), 
Aoo(a = 0.2), Aoo(a = 0.4). No arrows have been drawn for a = 0.1, because they 
would fall on top of the corresponding curves. 

the infection, and high values of the initial fraction of infectives asymptotically 
lead to the endemic state. For infectivities in (Ai, A2), there is an intermediate 
interval of reconnection rates where the system reaches the infection-free state. 
For sufficiently large values of f, though, the system returns to the bistable 
region. Finally, for infectivities smaller than Ai, the situation is the same as 
in the previously discussed models: increasing reconnection leads the system 
from the endemic state to the bistable region and, upon further increase of f, 
the system reaches the infection-free zone. 

Coming back to the question raised at the beginning of this section, about the 
possibility of reconnection having unwanted effects in the spread of the disease, 
the results shown in Fig. [7] suggest that the answer is not straightforward. As a 
function of the reconnection rate f, the stable equilibrium fraction of infectives 
can grow or decrease and, even for a given value of r and A, the sign of the 
derivative of nj with respect to f depends on a. As a partial characterization 
of this diversity of behavior, we analyze the sign of dnj/df for f — ► 0, i.e. 
the slope of the stable branches depicted in Fig. [7J for small values of the 
reconnection rate. Figure [8] shows the zones of the parameter plane (A, a) 
where the equilibrium fraction of infectives increases with the reconnection 
rate, for f — ► 0, and for two values of k. It can be shown that, nevertheless, 



14 




Fig. 7. Bifurcation diagram for the total fraction of agents carrying the infection, 
ni + ha, for a disease with an asymptomatic stage of duration a/7, for k = 3, 
A = 1 and several values of a. Full and dashed lines represent stable and unstable 
equilibrium solutions, respectively. Dots are simulation results for a system of 5000 
agents with 7 = 1. The arrows indicate the asymptotic value of ni for a = 0.4 and 
a = 0.8. 



the growth of nj within this regions is small if the infection level is itself low. 
In any case, in most of the parameter plane, the fraction of infectives exhibit 
the desirable decline of the infection level with reconnection. 



It is tempting to apply these results to the propagation of real diseases. The 
parameter a can be associated wit h the number of in fections that occur be- 
fore the onset of clinical symptoms (IFraser et all 120041 ). Using estimates of this 
number, we can have an idea of how contact switching mig ht work for these in 



fectio ns. For gonorrhea, a rough estimate, using data from (IKretzschmar et al. 
19971 ). gives a ~ 0.2 and A ~ 4. Being a sexually transmitted disease, k = 



is a reasonable assumption. Using these values, we see from the upper panel 
of Fig. [8] that gonorrhea would be in the zone whe re contact switchin g is not 
effective. For influenza, the estimate gives a ~ 0.4 flFraser et all I2004T ). omce, 
in contrast with gonorrhea, influenza is an airborne disease, the number of 
contacts that may be infected is large. Using k = 10, the lower panel of Fig. [H] 
shows that, almost for any infectivity, contact switching might be an effective 
strategy to curb the epidemic outbreak. Needless to say, any application of 
these ideas should be based both on more realistic models, taking into ac- 
count specific features of the disease under study, as well as on more sound 
estimates of the parameters involved. 
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Fig. 8. Phase diagrams showing, for two connectivities, the regions (shaded) in 
parameter space (A, a) where the derivative of the fraction of agents that have 
contracted the disease with respect to f, at f = 0, is positive. 

5 Discussion and conclusion 



We have studied the effects of contact switching as a strategy to prevent 
the spread of an epidemic disease. Our analysis confir ms, and extends in 



several directions, the results obta ined in previous works (j Gross et all 12006 



Zanette and Risau-Gusmanl . 120081 ). In particular, we have here relaxed some 
important constraints, for inst ance, the assumpti on that every agent knows 



the state of every other agent (IGross et al.l . 120061). and the assumption that 



the m ean number of contacts per agent is large (jZanette and Risau-Gusmanl . 
20081 ) . In the basic model that we analyzed first, susceptible agents with infec- 
tive neighbors are allowed to reconnect these links at rate r. For small values 
of this reconnection parameter, there is a continuous transition -through a 
transcritical bifurcation- from an infection-free state to an endemic state, as 
the infectivity grows. When r is larger than a certain critical value, on the 
other hand, the transition becomes discontinuous and gives rise to a bistabil- 
ity region where the infection can become asymptotically extinct or endemic. 
In this situation, the asymptotic state is determined by the initial fraction 
of infectives. Further increase of the infectivity, however, makes the endemic 
state globally stable, and the infection-free state becomes unstable. 



One crucial feature of our basic model is that only susceptibles are allowed to 
reconnect. Even though this constraint might seem reasonable, the voluntary 
nature of contact switching suggests relaxing it. We have analyzed a model 
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where both types of agents can switch contacts, with a new parameter p giving 
the probability that it is the infective agent who breaks and reconnects the link 
with a susceptible neighbor. Remarkably, we find that, for certain values of the 
parameters that characterize the disease, reconnection can completely suppress 
the infection for intermediate values of p, but not for larger or smaller values. 
Furthermore, even when the disease is not completely suppressed, equilibrium 
infection levels can be lower when more infectives are allowed to reconnect. 
Even though, at first sight, this might seem counterintuitive, it can be readily 
understood, at least for high infection levels: if the fraction of infectives is large, 
susceptibles are likely to reconnect their broken links to other infectives, thus 
not altering their risk of getting infected. On the other hand, when an infective 
agent reconnects, which with high probability occurs to other infectives, the 
number of "dangerous" links decreases. 



It has been shown that the outcome of many strategies to control epidemic 
outbreaks depends not only on the reproductive number of the disease, but 
also on the number of inf e ction s caused by an agent before the onset of clinical 
symptoms (IFraser et all 12004] ) . In the case of contact switching, it is even 
possible that the presence of an asymptomatic phase leads to an acceleration 
of the spread of the disease. In Section H] we analyzed an extension of the 
basic model, introducing a parameter a which gives, on average, the fraction 
of time that the agent spends in an asymptomatic phase during the infection. 
In this phase, the agent behaves as a susceptible individual, but can transmit 
the infection through contagion. We find that, in this case, reconnection does 
lead to higher levels of infection, but this effect is restricted to some limited 
regions of parameter space. Furthermore, within these regions, the increase 
of this infection level due to reconnection is rather small, unless its value is 
already very high in the absence of reconnection. 



It should be clear that any significant assessment of the benefits and risks of 
contact switching must use sound estimates of all the variables involved, ap- 
plying them to realistic models of the analyzed population. Important features 
that may be added to the above models are, for instance, the population divi- 
sion into females and males -who often exhibit rather different epidemiological 
responses- and age groups, as well as more complex degree distributions in 
the initial network. Such generalizations, however, will hardly be amenable to 
analytical treatment. In our basic model and its extensions, on the other hand, 
most epidemiological features can be studied analytically. It is only after these 
have been fully understood that more complex models should be tackled. 
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A Appendix 



Models of epi demic spreading can be s tudied as stochastic contact processes 
on networks (ILevin and Durrettl . Il996l ). We illustrate the procedure with the 
basic SIS model presented in Section [2j For each site x of the network we 
introduce the probability P t (A x ) that the agent at x is in state A (either S or 
I) at time t. The evolution equation for this probability is: 



P t {I x ) = -P t (S x 



y^x 



Iy) 



(A.l) 



where 7 is the rate at which an infective agent becomes susceptible, and A is 
the rate at which infective agents infect each of their susceptible neighbors. 
Here, Pt(S x ~ I y ) denotes the probability that the agents in sites x and y are 
susceptible and infective, respectively, and that they are connected. Generally, 
we also introduce Pt{A x ~ B y ) as the probability that an agent in state A at 
site x and an agent in state B at site y are connected, and Pt(A x 00 B y ) as the 
probability that the same agents are not connected. The evolution equations 
for these two-sites probabilities are 



P t (S X ~ Iy)= 1 P t (l X ~ /„) - 1 P t (S X ~Iy)+\J2 P t(Sx ~Sy~I z 

z^x,y 

x J2 p t(i y ~ s x ~ I z ) - XP t (S x ~ /„) + 

z^x,y 

P t (ly >*S X ~ I z ) _ P t {ly ~S x OO Q z ) 



r 



z+x , y Ev^Z Pt(Sy <* O v ) y J2v^X Pt(Sy *> O z ) 

P(S X ~ Sy)= 1 (P t (I x ~ Sy) + P t (S x ~ Iy)) - 

A( Pt(S x ~S y ~I z )+ Pt{S y ~ & ~ I z ) + 



v P t (S x ooSy- I z ) v P t (S, 00 S x ~ 2 



-PtCSi, "° 0«) E^x PsO^ ^ O z 

P t (l X ~ J„) = -2 7 P t (J, ~ J„) + A(P(S, ~ Iy) + (Sy ~ I X )) + 

A X! (^t(4 ~ S y ~ J Z ) + ~ ^ ~ J,)) (A.2) 

P(5 X . 06 Iy)=jP t (l X * I y ) - 1 P t (S x I y ) + 

A £ P t (^ OO Sy ~ I,) - A £ P(J, * ^ ~ J 2 ) - 
P t {Iy~S x *I z ) _ P t (4^^~0,) 



r 

z 



z ^ y J2v^z Pt(S y 00 O v ) J2v^x Pt(Sy Og) 

P t (l X ^Iy)=-2 1 P t (l x ^Iy) + 
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A £ (P t (l X *Sy~ I Z ) + P t {ly *S X ~ I Z )) 

where r is the rate at which a susceptible-infective link is broken, and the sus- 
ceptible agent reconnected to a random agent. The remaining two-site prob- 
ability Pt{S x oo S y ) is obtained from the fact that the sum of all two-sites 
probabilities equals unity. The interpretation of the different summands in 
each equation is rather straightforward. In the first equation, for instance, 
the first two terms account for the probability that a connected susceptible- 
infective pair is created from an infective-infective pair, or destroyed, by the 
spontaneous cure of an infective agent. The third and fourth terms account for 
the probability that infection is acquired from, or transmitted to, a neighbor 
outside the considered pair. The fifth term accounts for the probability that 
infection is transmitted from the infective agent to the susceptible agent in 
the same pair. The last two terms account for the probability that the pair is 
created or destroyed by the reconnection process. In these terms, the denomi- 
nators take into account the fact that, when the susceptible agent reconnects 
a link, the new contact must exclude the previous neighbors of the agent in 
question. 



Within the pair approximation (ILevin and Durrettl . 119961 ). the three-site prob 



abilities in Eqs. (1A.2I) can be approximately given in terms of two-site proba- 
bilities, thus closing the equation system. The basic idea of this approximation 
is that, in a chain of three of connected sites, x ~ y ~ z, the mutual interac- 
tions between the central site y of and each of the other two neighbors, x and 
z, can be decoupled: 



.P t {A x ~ B y )P t {B y ~ C z 
Pt(B y 



The constant K is determined either from analytical considerations, which 
give K — (k — l)/k, or from heuristic arguments -for example, a better fit of 
simulations, as done here in the main text. 

We assume that all the probabilities are homogeneous over the population, i. e. 
that they do not depend on the site, which is true if the initial probabilities are 
also homogeneous. We also assume that the epidemiological states at sites that 
are not connected can be taken as independent, i. e. Pt{A ^ B) « P t (A)P t (B) 
and Pt{A ^ B ~ B) w P t (A)P t (B ~ C) (using also the pair approxima- 
tion). A further approximation that simplifies the equations considerably is 
to assume that the number of neighbors of each site is much smaller than the 
total number of agents so that, when an agents reconnects, the probability 
that another specific agent becomes the new contact is simply (N — 1) _1 , thus 
simplifying the denominator of the reconnection terms. 
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To write the equations in terms of the fractions of agents, ns and nj, and 
links, msi , rriss and rrijj we take into account the relations 



n s (t)=P t (S), 
m ss (t) = ?^±P t {S~S), 



m S i(t) 



k 
2(N 



(A.4) 



r 



k 



Pt(S~I), 



and nj = 1 — n s , m n = 1 — m ss — mjj. Using these equivalences, and the 
approximations mentioned above, Eqs. ([!]) are obtained from system flA.21) . 
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